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model for the planar graphene sheet. With the model, the low frequency behavior of the dispersion 
of the acoustic modes as well as the flexure mode can be precisely calculated. Based upon a com- 
prehensive chiral symmetry analysis, the calculated mode frequencies (including all the Raman and 
infrared active modes), velocities of acoustic modes and the polarization vectors are systematically 
fitted in terms of the chiral angle and radius, where the restrictions of various symmetry operations 
of the SWCNT are fulfilled. 
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I. INTRODUCTION 

Since the discoveryJ. of carbon nanotubes, it stimulated 
extensive studies on the lattice dynamics of the single- 
wall carbon nanotubes (SWCNT) ,Vi 4 i 5 i 6 The SWCNT 
can be viewed as rolling up a two-dimensional graphitic 
lattice into a seamless cylinder^ Among various stud- 
ies, the zone-folding method was firstly applied to study 
the lattice dynamics of the SWCNT, which is migrated 
from that for the planar graphite sheet. It successfully 
accounts the general features of the phonon spectrum 
in SWCNT. However, unlike the electronic problems, on 
which this method is workable except for very narrow 
SWCNT, the results for the lattice vibration by the zone- 
folding were not very satisfactory, 2 particularly, it can 
not provide correct numbers of Raman and infrared (IR) 
modes of SWCNT£ 

In recent a few years, the theoretical investigation on 
lattice dynamics of SWCNT has achieved many success- 
ful and interesting results. Using the rotational and he- 
lical symmetries, Popov et al&HL decomposed for the 
first time the lattice dynamics equation of the SWCNT 
into a six-dimensional eigenvalue problem. Analytical ex- 
pression for the velocities of the acoustic modes are de- 
rived using Born's perturbation technique.— The trans- 
lational and rotational sum rules expressed in terms of 
the force constant matrijiii of the helical tube are also 
addressed in Ref.llOt Damnjanovic et al. applied the line 
group analysis to the lattice dynamics of SWCNT 
and also obtained the generalized Bloch form by modi- 
fied group projectors techniquei 14 i 15 > 16 They also calcu- 
lated the phonon spectra in SWCNT with a modified 
force constant modeiiLiSiiS in which certain kinematic 
constraints were imposed to achieve the twisting mode. 
The authors tried also to fit their calculated frequencies 
of the Raman and IR modes as functions of tube radius 
and chiral angel. In a series of papers Mahan and his 
collaborators^}^ worked on the lattice dynamics of the 



SWCNT with various vibrational energy terms such as 
spring energy, in-plane or out-of-plane bond bending po- 
tential energies successively, and calculated the phonon 
spectrum for the achiral as well as chiral 2 ^ tubes. They 
realized that, to keep the twisting mode to be a zero 
mode in the long wave length limit, all the pieces of vi- 
bration energies under their consideration violating the 
global rotational invariance around axial axis should be 
ignored. The authors also interpreted the existence con- 
dition of the flexure mode, which is an inherent character 
for rods and plates. This mode also shows up in Popov 
et. aVs calculation but is absent in Damnjanovic et a/.'s 
results. 

However, as pointed out in Ref. l2ll the commonly used 
expression for bond bending potential, even in addition 
with the kinematic constraints as in 

Refs|l2|and|H does 
not satisfy the rotational invariance condition, where the 
missing of the flexure mode is a kind of evidence. At the 
same time, the discussion on the rotational symmetry 
requirements to the potential terms in connection with 
the flexure mode is restricted only to those around the 
tube axis. The force modeliS taken from that for the 
graphite^ 3 , obeys the translational and rotational sum 
rules of the planar sheet but does not meet precisely the 
symmetry requirements of cylindrical warped hexagonal 
lattice sheet. 

Actually, the rotational invariance symmetry, i.e., the 
potential energy must keep unchanged when a system is 
rigidly rotated, should play a nontrivial role for rod-like 
cylindrical lattice sheet. On the other hand, the intrinsic 
symmetry of the tube lattice structure would also lead 
to specific restriction on dynamic matrix. For a com- 
plete account of the symmetries of the lattice dynamics 
for the SWCNT, both aspects should be considered in a 
consistent way. This constitutes the basic scope of this 
paper. 

As a continuation of the previously mentioned series 
of works, in this paper, we show the vibration potential, 
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i.e. force model, up in a detailed expression with the 
rigid translational and rotational invariance being pre- 
cisely guaranteed, in which the curvature effect on bond 
lengths and bond angles is carefully considered. As its di- 
rect consequences, the flexure mode as well as the acous- 
tic torsion mode can be correctly treated. Moreover we 
calculate the phonon spectrum with the proposed force 
model and carry out a full chiral symmetry analysis. In- 
stead of previous case-to-case empirical forms, 1 - the fre- 
quencies of all Raman and IR active modes, the polariza- 
tion characteristics of the modes, as well as the velocities 
of all four acoustic modes can be systematically fitted in 
terms of tube radius and chiral angle with the symmetry 
considerations. 

As expected, our calculated results show that the cur- 
vature effect on bond lengths and bond angles would 
be no more negligible only for nanotubes of smaller ra- 
dius. Meanwhile the nontrivial contribution for the 
dependence of varies physical quantities is also effec- 
tive only for tubes with enough small radius. For thin 
tubes with the diameter range such as 4 ~ 10 A, the 
warped graphene sheet model, on which our discussions 
are based, should be fairly applicable to the SWCNT. 
For tubes with diameter smaller than 4 A approximately, 
the warped graphene sheet model might have to be 
improved by further optimizing its equilibrium lattice 
configuration)^i2^2SiSi and the part of chiral symme- 
try analysis in this paper were not be valid any more. 
Such estimation is entirely consistent with the existing 
literatures i& 

The present paper is organized as following. In Sec. II, 
after a brief description on notations, we present the de- 
tailed expression of the vibrational potential. The rigid 
rotational invariance and the comparisons with previ- 
ous works are discussed. Subsection III A is devoted to 
a generic chiral symmetry analysis for various physical 
quantities, while the main results and relevant discus- 
sions on phonon spectrum calculations are presented in 
Subsec. Ill B. The paper ends with a brief summary in 
Sec. IV. 



II. MODEL 

A. Notations 

The z axis in this paper is set along the tube axis and 
the x axis comes across the middle point of a C-C bond in 
Fig. It is understood in this subsection that notations 
for cylindrically warped tubes are often introduced with 
the aid of those defined on the planar sheet. 

As is well known, SWCNT can be notated by a chiral 
vector R — nidi + 712(12 or equivalently, the radius r 
and chiral angle 9. Where di and 0*2 being the primitive 
lattice vectors in graphite lattices^ |ai| = |a2 1 = 2.46 A. 




FIG. 1: The sketch of carbon atoms on a cylindric surface. 
They have the relations as 



r = \ n~( + nifi2 + nz, = atan . 1 

2?r V 1 2 ' 2ni +n 2 

We use the helical and rotational symmetry descrip- 
tions for the SWCNT throughout this paper. Introduce 
a real lattice vector H = pidi -\-p2d2 describing the screw 
operation^ with nip2 — n 2 pi = N (N is the greatest com- 
mon divisor of ni and 712). The rotational angle of the 
screw operation is a = (H ■ R)/\R\ 2 (in the unit of 27r). 
A unit cell in the {n\, 77,2) tubule can be notated by two 

equivalent ways, either based on {H, or on (01,02) 
as 

-* R 

r m j = mH + l — , or 7^,^=0101+0202 (2) 
with the relations 

m = (77,102 - ?7 2 gi)/A, I = qipi - g 2 pi . (3) 

In the following, we use (H, J|) as the basic vectors, and 

bu and are their corresponding reciprocal unit vectors. 
Any wave vector in the reciprocal space can be written 
as 

K=^ H + ^b R . (4) 

Z7T iV 

In the first Brillouin zone, n G (— \, ^] and n is an integer 
in(-f,f]. 

B. Potentials 

It had been proposed^ that there are five distinctive 
terms for the vibrational potential of the graphene sheet 
satisfying the rigid translational and rotational invari- 
ance of the planar sheet. We herein present a detailed ex- 
pression of the vibrational potential for the SWCNT with 
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the curvature effect being carefully in-cooperated, which 
is essentially the generalization of that for the graphcnc 
sheet. It satisfies precisely the requirements of the rigid 
translational and rotational invariance and realizes the 
corresponding general symmetry sum rules in Ref. [ToL 
i.e., the potential energy must keep vanish term by term 
when the tube is rigidly translated or rotated around any 
axis. Introduce f as the equilibrium position of atom i 
and Ui as its displacement vector, fj = fj — f is the 
vector from atom i to j in the nanotube while the modu- 
lus represents the length of C-C bond between atoms 
i and j. f is determined following the geometry of a 
warped graphene sheet so that the three tridental bonds 
nearest-neighbored with the atom i as well as the angles 
between any of the two bonds are not equal to each other, 
especially for thin tubes. 

The five potential terms for the SWCNT are expressed 
in the following. 

(1) Vi is the potential of the spring force between the 
nearest- neighbor atom pair, 



k. 



2 ^ 

4=2 



[(ui - u\) 



If J 



(5) 



where kt is the first-order force constant and e, ■ = t^t- 

1 11 \ru\ 

We'd like to point out that the components of the dis- 
placement vectors perpendicular to e < li which violates the 
rigid rotational invariance are forbidden. 

(2) V s i is also the potential of the spring force but 
between the next nearest-neighbored atoms illustrated 
as (1, 5... 10) in Fig.[U 



(6) 



i=5 



with k s i the second-order force constant. 

(3) The potential energy for the in-surface bond bend- 
ing Vbb is actually a term associated with bond angle 
variations. Three atoms have to be considered together, 
and 



to an arbitrary axis can be kept only when the differ- 
ences among bond lengths and bond angles be carefully 
accounted. 

(4) The potential of the out-of-surface bond bending 
V rc is the energy between atom i and the three nearest- 
neighbor atoms ji with four atoms simultaneously con- 
tained, 



V r , 



■[(3i 



eH 2 , 



E 



(8) 



(9) 



Where i takes 1 or 2 with ji running over the three near- 
est neighbors of atom i. This potential has the physical 
intuition as that responsible for the radial optical mode in 
tubes and is trying to keep the four atoms on the cylinder 
surface. Herein we introduce an unit vector e£ c defined in 
Eq. @, which is rather close but not equal to the radial 
unit vector of site 1 shown in Fig. 2. When the radius 
is large enough, the e^ c would be close enough to the e[. 
However, we stress that the potential term V rc with e^ c 
keeps the rigid rotational invariance while that with e^ c 
substituted by e£ would break the rotational symmetry. 



raduis [angstrom] 



V BB 



iBB 



E E 



(4 



+- 



Uji — Ui 



'ij' 



(4k " CWfyijAj'.) 



<BB 



(7) 



3* fi 



If i is 1 or 2, take the sites 2,3,4 or 1,5,6 respec- 
tively. In Eq. (J7|, djiiji stands for the equilibrium angle 



between the bonds f 



and fiAi , while 6'. . ., 



for the corre- 



sponding angle in vibration. The expression in terms of 
u,j i — Ui and Uj>. — Ui is an honest representation to that 
of cos#', where the rigid rotational invariance referred 



FIG. 2: The projection of ei c on e[ for tubes (2n, n) with 
n £ [1, 15]. It shows that e^ c only deviates about 2% from e[ 
even in the small radius (2, 1) tube. 

(5) The twist potential energy for bond f\h is 
V tw = ~y E [("< - - - %')) ' elfc] 2 » ( 10 ) 

<i,j> 

where e[ k is the unit vector along the radial direction 
of the middle point of fifc, < i,j > represents a pair of 
atoms nearest-neighbored with atom 1 while k the third 
of its nearest neighbors. Pair < i',j' > is the image of 
< i,j > referring to a C2 rotation around the axis in 
e[ k . It has the intuition that responsible for modes with 
twisted vibrations. 
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C. Some details about the model 

Obviously, all the five potential energy terms respect 
the translational invariance. When Uj = Uj, u.i — uj = 
leads to 

Vl = Vsi = V BB = V rc = V tw = . 

For the rotational invariance, we have to consider the 
five potentials term by term separately^i^ When the 
tube rotates rigidly around an arbitrary axis for a small 
angle Suj with its direction along the axis ||^y and its 

magnitude \5u)\ being the totation angle, each lattice site 
acquires a displacement Ui — Suj x fi, 

Ui — Uj = Suj x (fi — fj) = Suj x fji . (11) 

Substituting Eq. (|11|) into the first two potential terms 
© and ©, it is straightforward to have (uj — Ui) ■ eL = 
r„(fc;x e i.)-e1 3 =0. Then 

Vi = V s i = 0. 

Substituting Eq. I|ll|) into the third potential term , a 
typical representative term in summation becomes 

v BB ~ • (4 >< 4 + 4 x 4 2 )] 2 = o . 

In which a fact has been used that r%j in the denominate 
is cancelled by that in the numerator when Eq. I|ll|) is 
applied. Moreover, for each typical representative term 
in potentials (jHJ and l|10(l . we have similarly, 

V„ ~^[Sl3x (r 12 + f 13 + r 14 ) • e[ c ] 2 = , 
V tw ~^[5ujx (r 43 - r 56 ) • e^] 2 = . 

As an example for further clarification, one might cal- 
culate the phonon spectrum with all the bond lengths and 
bond angles assumed to be equal to that of the graphite. 
It results that the twisting mode (TW) at (k, n) = (0, 0) 
is no longer a zero mode and there is a finite gap. For the 
SWCNT (5,2), the error introduced by the equal-bond- 
length creates a gap of the order 0.5 cm -1 as shown in 
Fig. |3J where the radius r = 2.45 A, the three C-C bonds 
are about —1.3%, —0.3% and 0.0% shorter than that of 
graphite 1.42 A respectively. Although it is a minute 
number and entirely negligible in practice, it is of quali- 
tative significance. 

As a distinguished feature of the nanotubes that the 
two degenerate transversal acoustic (TA) modes shown 
up at (k,u) = ±(a, 1) are flexure modes. Instead of the 
conventional linear behaviors, the low frequency limits 
of their dispersions are parabolic as uj 2 = /3 2 (k =F a) 4 . 
We address that the rigid rotational invariance around z 
axis itself is not a sufficient condition for the existence 
of the flexure modes cf. Ref. |2lJ As a counter exam- 
ple, we introduce a term V Ti = ^f-[(u2 ~ ui) ■ e^] 2 with 
^12 — ^12 x ^12 : which is the elastic energy due to the 
tangential spring force between atoms 1 and 2. It is easy 




FIG. 3: The effect of bond lengths on TW mode in tube (5, 
2). (a) All bonds are assumed to be the same. The frequency 
of TW mode is nonzero, (b) The differences between bonds 
are considered. The frequency of TW mode is presicely zero. 

to verify that V T i = when the tube rotates around the z 
axis. But there is no flexure mode when V T i is taken into 
account^ This is because that V T i would be no more zero 
when the tube rotates around any axis perpendicular to 
the z axis. 

III. RESULTS AND DISCUSSIONS 

A. Further symmetry considerations 

We may establish a mapping from the space of chi- 
ral vectors on the planar sheet to that of the nanotube 
structure in a fixed way of wrapping 29 

/ : {chiral vector set} i— > {nanotube set} . 

However this is not a one-to-one mapping. As mentioned 
m Ref.lH a given SWCNT can be equivalently composed 
by three chiral vectors Rq, R2 and i?4 with chiral angles 
as 9, 8 + 27r/3 and 9 + 4ir/3 respectively within the cor- 
responding graphite sheet. Due to the Cq v symmetry of 
the hexagonal graphene lattice sheet, there is another set 
of vectors Ri, R3 and R$ with chiral angles as 9 + tt/3, 
9 + it and 9 + 5ir/3 respectively on the graphite sheet. 
These three chiral vectors correspond again to the same 
SWCNT which is actually the nanotube by rotating the 
tube formed by Rq upside down, i.e. a Cix operation. 
The net result of the operation is an exchange of A and 
B carbon atoms in unit cells with the sign of its chiral 
index v = mod{ni — ?i2, 3} also changed. There is one an- 
other operation a xz connecting a pair of SWCNT which 
is the mirror reflection onto each other with respect to xz 
plane. Correspondingly Rq in the sheet is changed into 
its countpart R with 9^—9 but v kept unchanged. 

As one of the direct consequences of the above obser- 
vations, any physical quantity Q"(r, 9) of the SWCNT is 
a periodical function of the chiral angle 9 with period ^ , 
i.e., Q"(r,9) can be expanded asS 
00 

Q"(r, 9) = J2 a n cos(3n8) + b v n sin(3n0) . (12) 

Furthermore, under the operation 9^9 + 7r/3, any 
scalars S would keep unchanged, For normal vectors, the 
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radial components have the same behavior, but the az- 
imuthal and axial components change signs— That is 



S"(0+§) =S- V (6) , 
V J {0+l) = C 2x v- v {9) 



(13) 



with C 2x e r = e r , C 2x e$ = — e^, C 2x e z = -e z . Here e r , 
and e z are unit vectors oriented towards radial, az- 
imuthal and axial directions upon the cylindrical surface 
respectively. 

In addition, operation a xz gives that any scalars or 
radial and axial components of vectors are even functions 
of 9, while azimuthal components are odd functions of 
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S u (-6) = S v {0) , 
v»{-6)=a xz v»(0) 



(14) 



with a xz e r = e r , a xz e,p = -e^, & xz e z = e z . The detailed 
expressions reduced from Eq. I|12fl due to the symmetry 
restrictions Eqs (| 1 31) and 114fl are shown in the appendix. 



B. Frequencies, acoustic velocities, and 
eigenvectors 





We calculated the phonon spectrum in the (k, n) rep- 
resentation with the vibrational potential Eqs (|5|)- (|10fl . 
By tuning the calculated results to the experimental 
data— for the Raman modes of (10, 10) tube (see Ta- 
ble UJ , we fit the corresponding force constants as ki — 
364.0 Nm" 1 , k sl = 62.0 Nm" 1 , k BB = 1-07 x 10" 11 erg, 
k rc = 14.8 Nm" 1 , and k tw = 6.24 Nm" 1 . We successfully 
identified 4 zero modes, 14 Raman active modes and 6 
IR active modes for the chiral SWCNT£ Two of the six 
modes at (k, n) = (0, 0) are zero modes. They are lon- 
gitudinal acoustic (LA) and TW modes both belonging 
to o^o representation (Reps). Another mode belonging 
also to o^o Reps is the IR active e r optical (OP) mode 
with A and B atoms oscillating out of surface in radial di- 
rection. The other three nonzero modes at («, n) = (0, 0) 
belonging to qAa Reps are Raman active, they are e r 
acoustic (AC), eg OP and e z OP. There are one flexure 
and five nonzero modes at (k, n) = (a, 1). The latter are 
both Raman and IR active assigned to a E\ Reps. At 
(k, n) — (2a, 2), six nonzero modes are all Raman active 
belonging to the same 2 aE 2 Reps. 

We also calculated the sound velocities of LA and TW 
modes, the flexure parabolic coefficient for the TA mode, 
and polarizations of nonzero modes at (k, n) = (0, 0) for 
varieties of tubes as the functions of r and 9. The modes 
at (k, n) — (a, 1) as well as (2a, 2) are doubly degen- 
erated with their one-to-one degenerate counterparts at 
(— a,~ 1) and (—2a, —2) respectively. In Fig. 0| we pro- 
vide the three-dimensional plots for some of these quan- 
tities as functions of r and 9. 

Our calculated results arc further fitted following the 
parameterization shown in Eq. 1)12(1. The correspond- 
ing expressions for the frequencies of above 15 Ra- 
man and IR active modes, the velocities and flexure 




FIG. 4: As functions of r and 8, we show in (a) the frequency 
of e z OP mode at (ft, n) = (0,0); (b) C T w for small radius 
tubes; (c) the z component of the eigenvector of e r AC mode. 



parabolic, and the polarization vectors of nonzero modes 
at («, n) = (0, 0) are listed in Tables GH HTU and M\ re- 
spectively. In this paper particular emphasis is drown 
to tubes with smaller radius. So we take the parameter 
regime for frequencies and velocities in r S [4.0, 10.0] A 
and 9 £ [— fjf], while that for polarization vectors in 
r G [2.0, 10.0] A and 9 £ [—§,§]• We keep the relative 
errors in fitting less than 5 x 10~ 4 . The above parameter 
regimes are such chosen only for tubes of smaller radius, 
in which the contributions of 9 dependence are notable 
comparing to those of r dependence. The velocity of 
the TW mode (in Table ITTTf) as well as the polarization 
vectors of modes with o^o~ Reps (in Table Hv|) are typi- 
cal examples shown up an evident 9 dependence. In fact, 
the corresponding calculations have also been carried out 
in the parameter regime r € [3.13,58.7] A. The results 
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TABLE I: Comparison between the calculated results (following our model) and the experimental values for several mode 
frequencies (in the unit of cm" 1 ) of SWCNT (10,10). 



Reps 






aE 1 


a E 1 


2aE^ 


2a E 2 


2a E 2 


2a E 2 


theory 


167 


1588 


105 


1588 


21 


367 


873 


1584 


experiment 30 


186 


1593 


118 


1567 


/ 


377 


/ 


1606 



TABLE II: Frequencies (in the unit of cm 1 ) of 15 Raman and IR active modes as functions of r (in A) and 6. 



Reps 


Mode/Order 


ui{6) 








Mr) 




oA+ R 


e r AC, 1 


Mr) 








Mr) 


iioo.oo loy.Do 




eg OP, 2 


Mr) 


+ Mr) 


cos 66 




Mr) 


- 1594.00 26%98, fl ( r ) - ^ 




e z OP, 3 


Mr) 


+ Mr) 


cos 66 + f 2 (r) 


cos 126 


Mr) 


1 t- r\ a r\r\ Q1 81 r / \ 1 S 68 r 1 \ fi8 

= 1594.00 - ^^f 1 , = — ^r-, M r ) = -^73^ 




e r OP, 4 


Mr) 


+ Mr) 


cos 66 




Mr) 


^ n a r~.t QQft 99 1117 ^0 c / \ Q Ifi 

= 864.81 + 22^22 _ / x (r) = ^ 


aEi R Ir 


1 


Mr) 


+ Mr) 


cos 66 




Mr) 


(±u.±o 1 410 . U( / „ \ 1 . yy oi.oo 




2 


Mr) 


+ Mr) 


cos 66 




Mr) 


_ 1603.51 746.51 y ^ _ 115.54 




3 


Mr) 


+ Mr) 


cos 66 




Mr) 


^ o j , 860 Of) 1 758 70 r / \ 1163 206 52 

— 864.84 H p— — ^ — , /i(rj - -p- t .4 




v = ±1, 4 


Mr) 


±/i(r) 


cos 96 




Mr) 


" 1594.13 /i(r)-2L|2 




v = 0, 4 


Mr) 


+ Mr) 


cos 66 + f 2 (r) 


cos 126 


Mr) 


= 1594.14 - 21^48, / x ( r ) = L83 _ 12^)3 

p r \ 9 7n n fin 

Mr) = + ^ 




5 


Mr) 


+ Mr) 


cos 66 




Mr) 


- 1593.97 ™, ^(r)- 


2a E2 R 


1 


Mr) 


+ Mr) 


cos 66 + Mr) 


cos 126 


Mr) 


_ 959.33 736.60 | 779.59 

/i(r) = ^ + ^,/ 2 (r) = -2^ + 2^ 




2 


Mr) 


+ Mr) 


cos 66 




Mr) 


1420.21 , 54.52 1246.29 r { „\ 204.34 
1 r 3 r 5 i Jl( r ) r 'i 




3 


Mr) 


+ Mr) 


cos 66 




Mr) 


2535.48 2426.65 ( ( „\ 412.23 
r r '.i ! Ji( r ) r a 




4 


Mr) 


+ Mr) 


cos 66 




Mr) 


= 864.80 + ^1 471^81 + 12425JU 

Mr) = ^ ^ 




5 


Mr) 


+ Mr) 


cos 66 




Mr) 


- 1594.OO 269^19 + 22^21 
/r« = + 22|#i 




6 


Mr) 


+ Mr) 


cos 66 




Mr) 


- 1594.01 3»™1 2160,15 + Ml%26 



Mr) = H 



can be fitted very well by the same functions with minor 
change of coefficients, i.e., our fitting formulae could be 
applied to a much wider parameter regime. 

It is interesting to see that all the numerically fitting 
expressions meet the symmetry requirements of Eqs <(Kl|) 
and (|14fl . where the velocities and frequencies are of 
scalars while polarizations are of vectors. We notice that 
in Table [n] one of the a E\ modes manifests different pa- 
rameter dependence for different chiral index v. i.e. it 
has different expressions for v = ±1 and v = respec- 
tively. However, it is still consistent with the general 
expressions Eqs (I13|l and i|14fl . We further notice that, 
different from those physical quantities of normal vectors, 
the polarization vector can be measured up to a global 
phase factor as ±1. Therefore, the corresponding trans- 



formation properties with respect to 8 — > 9 + -| as well 
as 9 — > — 9 need to be generalized to 

v v J m \e + |) = \^\C 2x )C 2x v- v{m) {8) , 
v»; m \-9) = \(™)(a xz )a x yj m) (0), 

respectively, in which A 1 ™-* (6) is the phase factor taking 
value 1 or —1 corresponding to mode m such as e r AC, 
OP and e z OP and referring to the symmetry operation 
6 such as C 2x and a xz . In particular, we have 

X^ Ac ^C 2x )= X^ AC \a xz )= 1, 
\(^op) { c 2x) = A (eVOP)(^) = _i, 

A» OP )(C 2s ) = -A(«- OP )(a M ) = -l. 
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TABLE III: Sound velocities (in kms L ) of the TW and LA modes, and /3 (in 10 6 m 2 s 1 ) of the flexure mode as functions of 
r (in A) and 8. 

Velocity Velocity (8) /i(r) 

C TW / o (r) + /i(r)cos60 /„ = 13.5 - / x = ^ 

Cla /o(r) + /i(r)coe6fl /o = 21.0706 + « _ owao _ ^ = 002221 _ Mlfl 

P Mr) + /i(r)cos60 / = 1.3767r - 0.00142r 2 - 5.8 x 10 _5 r 3 , /i 



TABLE IV: Polarization vectors u = (u(A),u(B)) as functions of r (in A) and 8. Where u(A) and u(B) indicate the displacement 
vectors of atoms A and B in the (0, 0) unit cell respectively. For the three modes in this table, u r (B) = u r (A), u<t,{B) — —u r / > (A), 
u z {B) = -u z {A). 



(K,n) = 


(0,0) 




Vector(0) 


Mr) 






Ri (e r 


AC) 


U r (A) 


Mr) 


Mr) 


= 0.7071 

r z 








U4>(A) 


fi (r) sin 38 


Mr) 


0.0518 1 0.0468 
— ~r 5 — 








u z {A) 


fi (r) cos 38 


Mr) 


0.0517 1 0.0749 
\ r 2 




R2 {e<f> 


OP) 


u r (A) 


fi (r) sin 38 


Mr) 


0.0542 0.0455 
r r z 








u<t,(A) 


Mr) + Mr) cos 128 


Mr) 


= 0.7056 + Mow, /l(r) = 


0.0015 - 






u z {A) 


/i(r)sin 66» + / 2 (r)sinl26> 


Mr) 


- 0.0656 ™ )/2(r) - 


0.0048 ^1 


R 3 (e z 


OP) 


U r (A) 


fi (r) cos 38 


Mr) 


0.0447 0.0417 
r 








U<f,{A) 


Mr) sin 66» + / 2 (r)sinl26> 


Mr) 


= -0.0656 + / 2 (r) 


- 0.0048 + Mill 






u,(A) 


/o(r)+/i(r)coBl2fl 


Mr) 


= 0.7056 + 2^1£,/x(r) = 


0.0015 MOM 



It can be easily verified that the parameterized expres- 
sions for the modes e r AC, OP and e z OP listed in 
Table Hvl satisfy exactly Eq. I|15ll . Interestingly, the val- 
ues taken by A^(6) coincide precisely with the signa- 
tures of the major component of the corresponding modes 
with respect to the operations and a xz . The detailed 
expressions for polarization vectors are given in the ap- 
pendix. 

It is known that the planar graphene is the r — > 00 
limit of the SWCNT. We notice that six modes in Ta- 
ble[n]evolve to the two degenerate in-plane optical modes 
of the graphene with approximately the same frequency 
limit 1594.0. While three modes with / (r) -> 864.8 
approach the out-of-plane optical mode of the graphene 
sheet. Moreover shown in Table ITTll the sound velocities 
of the two zero modes, i.e. LA and TW modes belonging 
to qAq , have nonzero limits with different values. There- 
fore, it is expected that those two modes would approach 
two acoustic modes of the graphene sheet. We further 
notice that the frequency of the e r AC mode belonging 
to qAq (Table |nj| vanishes with its polarization vector 
perpendicular to the limiting sheet. We might interpret 
that this mode is a kind of precursor of the flexure mode 
of the graphitic sheet. However, at (n,n) = ±(a, 1), the 
dispersion of the TA branch is quadratic in k. We stress 
that the parameterization for the coefficient (3 cannot be 
extrapolated to r — > 00. This is prohibited by a kind 
of symmetry argument that the rod-like tube has two 



flexure modes with the cylindrical symmetry while the 
plate-like graphene sheet breaks the symmetry so as to 
have only one flexure mode. There is no way to cross 
continuously from the former to the latter. 

The sound velocities, slopes of the dispersions, of 
zero modes can be determined correctly only if we have 
enough data of the dispersions with high precision in the 
low frequency limit around [k, n) = (0,0). As one of 
the immediate consequences of our calculations, besides 
that for C LA , our calculations provide a credible value 
for C TW by taking the advantage that our model keeps 
exactly the translational and rotational invariance. With 
the velocities C LA and C TW we can extract the chirality- 
dependent expression for the Poisson ratio fj,i& Ignoring 
-j and higher order terms, we get 



fi = 0.21 



0.206 0.43 



cos 60 



(16) 



We may take the values of 9 as either ^ (armchair) or 
(zigzag) in the above Eq. I|16fl to extract the Poisson 
ratios for the achiral cases: 



/' : 
/' : 



0.21 
0.21 



0.636 
0.224 



for e=\ 

6 

for 6 = 



The corresponding results reported in Ref. 10 derived by 
the perturbation approach can be fitted by our expression 
in good agreement. 
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IV. CONCLUSIONS = va + cti cos 30 + va 2 cos 60 + . . . . (A2) 

In conclusion, we presented a detailed expression for Moreover, for polarization vectors, when X m (a xz ) = 

the vibrational potential of the SWCNT. On the one ~1> and A m (C 2K ) = -1, we have 
hand, they obey precisely the rigid rotational and trans- 

lational symmetries, in which the bond length differences v(as _ [ 1+( ~ 1)n 
as well as those of bond angles are properly considered 

following the geometry of cylindrically warped graphene , . , . 

sheet. On the other hand, they exhibit the chiral symme- ~~ 1 Sm + ^ 2 sm + . . . 

try inherited from the planar hexagonal lattice. With this v [(j) t 1- '? 1 '" 1 3 q 

model we calculated the phonon spectrum of the SWCNT ^ ' ^ " 

and parameterized the physical quantities in terms of _ ™~ _)_ vai cos 3^ _j_ fl2 cos qq 
tube diameter and chiral angle to make the chiral sym- 



v" r {6) = — s — ib n sin3n9 

n=0 



00 

metry to be explicit. Our calculated results are mainly v v {9) = 'S^ v ^ ^^ ^b„sin3re0 

illustrated and summarized in Fig. 4 and Tables [nj lllll 7 



n=0 



and llVI The physics of relevant quantities such as flexure _ v ^ g j n -^q _|_ j, 2 sm qq _j_ (A3) 
modes, some of the immediate consequences like Poisson 

ratio, and the symmetry properties of, for example, the when X ™(a xz ) = 1)&ad X ™{C 2x ) = -1, we have 
polarization vectors are discussed. Our paper provides a 

comprehensive and systematic understanding for the lat- 00 

tice dynamics of the SWCNT, and forms a kind of sound v"(6) = ^ — 5 — 'a„cos3n0 
basis for further improvements and developments. 



71=0 



vclq + a± cos 39 + va 2 cos 69 



APPENDIX A: THE EXPRESSION FOR v "(q\ _ J 1 ^T^ lfe siri3nfl 

PHYSICAL QUANTITY 



= vb\ sin 39 + b 2 sin 69 + . . . , 

With the symmetry restriction as Eqs i|13fl and Q14[l. 00 

i.e. \ {m \C 2x ) = \ {m X<r xz ) = 1 in Eq. dJ), the general v v z {6) = ^ ^^^W cos3n9 

expression Eq. ltT^|) would reduce to the following. n —o 

Since the frequency u is a scalar, we have = a + va\ cos 30 + a 2 cos 60 + . . . . (A4) 

w±(0) = a ± ai cos(30) + a 2 cos(60) ± a 3 cos(90) + . . . ; When A m (<7 xz ) = -1, and X m (C 2x ) = 1, we have 

(Al ) 

w°(0) =a a +a 2 cos(60) + . . . . 



The three components of a vector can be expanded as 



v»(9) = ^v [1 - L ^ r - ] b n Bin3n9 

71=0 

oc = vb\ sin 39 + b 2 sin 69 + . 

<(0) = v [1 ] a n cos3n9 ™ r i+(-i)« 1 



a + fai cos 39 + a 2 cos 60 + . . 



71=0 

vao + <Xi cos 39 + va 2 cos 60 



= ^ [i± ^ 1 &„sin3n0 ~ r i+(- iy , 

0W ^ <(0) = V i/[-VH 6nSin3n 



n=0 ' 

b\ sin 30 + ^62 sin 60 + . . . , 



v z( e ) = xy i±i ^ la " cos3r^6, 



n=0 

h sin 30 + vb 2 sin 60 + ... . (A5) 
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